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We propose a new and unique dark matter candidate: ~ 100 eV to 

~ 10 keV sterile neutrinos produced via lepton number-driven resonant MSW 
(Mikheyev-Smirnov-Wolfenstein) conversion of active neutrinos. The requi- 
site lepton number asymmetries in any of the active neutrino flavors range 
from 10 -3 to 10 _1 of the photon number - well within primordial nucleosyn- 
thesis bounds. The unique feature here is that the adiabaticity condition 
of the resonance strongly favors the production of lower energy sterile neutri- 
nos. The resulting non-thermal (cold) energy spectrum can cause these sterile 
neutrinos to revert to non-relativistic kinematics at an early epoch, so that 
free-streaming lengths at or below the dwarf galaxy scale are possible. There- 
fore, the main problem associated with light neutrino dark matter candidates 
can be circumvented in our model. 

PACS numbers: 04.90.Nn; 14.60.Lm; 97.60.Lf; 98.54.-h 
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Ample evidence shows that the dominant component of the matter content of the universe 
is non-baryonic Dark Matter contributing a fraction of the critical density Q m ~ 0.2 to 1 
In this letter we suggest a novel means to allow sterile neutrinos with masses ~ 100 eV to 
~ 10 keV to be the Dark Matter. We note that the existence of light sterile neutrinos may be 
an inevitable conclusion if the neutrino oscillation interpretations of the Super Kamiokande 
atmospheric neutrino data, solar neutrino data and the data from the Los Alamos Liquid 
Scintillator Neutrino Detector are correct |§. 

We envision production of our sterile neutrinos from active neutrinos via a net lepton 
number-driven MSW (Mikheyev-Smirnov-Wolfenstein |Q) resonant conversion process oc- 
curing near or during the Big Bang Nucleosynthesis (BBN) epoch. Only neutrinos which 
evolve through resonances adiabatically are efficiently converted. In turn, lower energy neu- 
trinos have larger ratios of resonance width to oscillation length and, hence, evolve more 
adiabatically. Therefore, the resonant process can produce a final energy distribution for 
the sterile neutrinos which is grossly non-thermal and skewed to very low energies. These 
neutrinos have energy distribution functions which are "cold" . This, coupled with their rel- 
atively large rest mass, implies that they go non-relativistic at quite early epochs, therefore 
essentially becoming Cold Dark Matter (CDM), albeit with a cut-off at or below the dwarf 
galaxy scale. Below this cut-off scale, free streaming of these neutrinos can effectively damp 
out density fluctuations. Since this cut-off scale is much lower than any previous light neu- 
trino dark matter models (e.g., Hot Dark Matter, HDM 0, Warm Dark Matter, WDM [§, 
and mixed Cold plus Hot CHDM |7j models), we deem these sterile neutrinos to be "Cool 
Dark Matter" (CoolDM). 

The contribution to the matter density of our universe from an active neutrino species 
(z/ e , or u T ) with mass m v is 

"~ 91.5/i 2 eV [ ) 

Here h = H /(100 km/sec/Mpc) where H is the Hubble constant. Therefore, active neutri- 
nos with masses in the O(10) eV range naturally serve as Dark Matter candidates. They are 



"hot" because their kinetic temperature is relatively high as a result of their light masses. It 
has been shown that models with HDM as the sole Dark Matter ingredient of the universe 
fail to form galaxies at an early enough epoch. This is due to the large free streaming of 
light neutrinos which damp out density fluctuations below the free streaming scale [H] 



This generic problem of the large free streaming length of light neutrinos is difficult to 
circumvent. One possibility is to have non-gaussian "seeds" in the primordial density fields 
||. Another is to have light neutrinos comprise only a small fraction of the dark matter 
(Q u /Q m < 0.2) [0,0. 

WDM models consisting of O(10 2 ) eV sterile neutrinos have a free-streaming length 
much shorter than that of HDM |10| . Because of this free-streaming, WDM models have a 



reduced amplitude of density fluctuations at the galaxy cluster scale, which fits observations 
of structures at this scale better than the COBE-normalized standard CDM (sCDM) model 
(with Q m = 1 and h = 0.5) |TD[. However, density fluctuations at smaller galaxy scales 
(~ lh^ 1 Mpc) in WDM models may be too small []TU[ to accomodate observations of Damped 
Lyman-a Systems at high redshifts, and observations of "lyman-break" galaxies which are 
highly-clustered galaxies at redshifts z ~ 3 Jll| . 



The neutrino Cool Dark Matter candidate proposed here is "colder" (i.e., with a much 
smaller free-streaming length) than all three classes of models mentioned above. In particu- 
lar, the cold non-thermal energy distribution of our candidate sterile neutrinos causes them 
to move much more slowly than WDM sterile neutrinos with similar masses. The small 
free streaming lengths in CoolDM models enable a circumvention of the generic problem 
associated with light neutrino dark matter candidates. 

The production mechanism of our sterile neutrino CoolDM candidate is very different 
from that of its closest kin - the sterile neutrino WDM candidate. In the model proposed 
by Dodelson and Widrow sterile neutrinos were produced by active-sterile neutrino 
oscillation at the BBN epoch. The oscillation was driven by finite-temperature matter effects 



T2|| . The end-result was a sterile neutrino population with an energy spectrum resembling 
that of an active neutrino species but with a overall suppression in normalization |J. The 
average energy of these sterile neutrinos is therefore the same as the active ones, ~ 3.151T, 
where T is the temperature of the active neutrino species. In our model, the sterile neutrinos 
are produced by a resonant active-sterile neutrino transformation driven by a pre-existing 
lepton number asymmetry. The resonant transformation is only adiabatic at the low energy 
portion of the neutrino energy spectrum. This feature creates a non-thermal sterile neutrino 
population, with only the low energy states populated. 

We quantify this production mechanism in the density matrix formalism |13|| . In this 
formalism, the v a <-> v s (a = e, // or r, and v s is the sterile neutrino) system is described 
by a four vector (P , P). The number densities of the neutrinos are n Ua = (Po + Pz)/2 and 
n Vs = (Po — P«)/2- The evolution of the mixing system in the BBN epoch satisfies the 
equation: 

P = V x P + P z - DP ± , (3) 

where Pj_ = P x ~k + P y y. The D term is the quantum damping term that acts to reduce the 
mixing system into flavor eigenstates and suppress the oscillation. Numerically D ~ G 2 F T 5 , 
where Gf is the Fermi constant. The vector V is the effective potential of the oscillation. 
Its components are 

V; = ^sin2£, V y = 0, V z = - 6 -^cos29 + V a L + V^ (4) 
where 5m 2 = m 1 Vs — m 2 ~ m l s i @ is the vacuum mixing angle, and E is neutrino energy. 



The matter-antimatter asymmetry contribution to the effective potential is [12 



« 0.35G F T 3 [L + 2L Ua + £ L Vg ] (5) 

where L represents the contributions from the baryonic asymmetry and electron-positron 
asymmetry, ~ 1CT 10 ; and L„„ is the asymmetry in the other active neutrino species For 
convenience, we will denote £ = 2L Ua + Yl a ^p^p- The contribution to V from a thermal 
neutrino background is Vj 10 2 G|PT 4 fjj. 



Similar quantities for the anti-neutrino v a <-> v s can be denned and they satisfy 

p = v x P + P z - DP ± . (6) 

where V x = V XJ V y = V y , V z = -5m 2 cos 26 /2E - V£ + Vj, and D » D. 

The sign of £ does not matter in the Dark Matter problem. For definitiveness we assume 
that it is positive, so that V£ is positive. Then the u a <-> u s oscillation system encounters a 
resonance (V z = 0, or equivalently, a maximally matter-enhanced mixing) as a result of the 
non-zero £ at a temperature 

T~ - 9 ( ( £) MeV, (7) 



,10 2 eV/ VOX 

where e = E/T. At T ~ T res , V^J is completely negligible. We have implicitly taken cos# ~ 1 
since the mixing angle must be small. The anti-neutrino v a «-> v s oscillation, on the other 
hand, has \V Z \ > \ — 8m 2 /2E\. It is therefore suppressed relative to a vacuum oscillation 
situation. The net effect is destruction of C, because v a is resonantly transformed into sterile 
neutrinos while u a is not. 

Quite similar to the putative matter-enhanced resonance transition in the solar neutrino 
problem, the adiabaticity condition of the v a — > v s resonant transformation at the resonant 
energy bin e res is 

-l 
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> I- (8) 



In the equation, \V X \ is the transformation rate, |V^(de res /dV^)| is the energy width of the 
resonance, and |de res /d£| is the the speed of movement of the resonance energy bin across 
the neutrino spectrum as v a neutrinos of different energies encounter resonance at different 
temperatures. (Strictly speaking, our expression for the rate and the energy width of the 
resonances are only valid when collisions are not important, i.e., D \V X \. However, the 
adiabaticity condition is nonetheless the same when D ^> \V X \ because the suppression in 
rate due to collisions is compensated by the increase in the energy width of the resonance.) 
From Eq. (§]) we have 



V, 



dCres 

dK 



~ e res sin 29. 



(9) 



From Eq. (I7|), we have 



dt 



AH 



d£/dt 
£ 



(10) 



where H pa 5.5T 2 /m p i is the Hubble expansion rate and m v \ pa 1.22 x 10 28 eV is the Planck 
mass. The adiabaticity condition is therefore satisfied if 



4 x 10 



1/2 , £ n 3/4 



e^/ 4 sin 2 20 



1 



1 - (d£/dt)/4H£ 



> 1. 



(11) 



,10 2 eV7 V0.1, 

Assuming adiabatic neutrino evolution through resonances, we have d£/dt = 
/(e res )de res /dt = /(e res )e res [4iJ — (dC/dt)/ C] where /(e) is the neutrino distribution function. 
We therefore derive £ = C imt — Jq icb /(e)de where superscript "init" indicates initial values. 
As a result, \(d£/dt)/£\ < H is always true unless \£\ <^ |£ mit |. This implies that the 
adiabaticity condition Eq. (|ll]) holds true for u a v s vacuum mixing that is not too small 
(sin 2 29 > 10~ 9 ), until most of C is destroyed by the resonant v a <-> v s transformation. The 
final £ is ~ after the resonant conversion process. In this limit, the total change of the v a v a 
asymmetry is AL^ pa £ mit /2. This change is entirely due to the v a to v s transformation, 
implying that the v s sea produced has a number density that is a fraction 



(12) 



of the number density of an active neutrino species. The sterile neutrino contribution to the 
matter density today is then 



a. 



( m v 



h 



(13) 



,91.5/i 2 eV/ V343eV/V°- 5 / V 0A J 

with a, (3 — e, fi, r. Here we have assumed that all neutrino asymmetries are their initial 
values and we have dropped the superscript "init" (L v will not change anyway). 

From Eq. (|4|) and ( p~3| ) we note that if the sign of 2L Va + J2/3^ a L U0 is negative, the anti- 
neutrino counterpart of v s will be resonantly produced from u a , and may therefore become 
the Dark Matter candidate if the counterpart of Eq. ( |l~3"l) is satisfied. 



As the transition is only adiabatic when £ is a significant fraction of its initial value, only 
the low energy sterile neutrinos (which encounter resonances first while £ is still relatively 
large, see Eq. J7|) are resonantly produced. Higher energy neutrinos go through the resonance 
later, when most of L is damped. Evolution through resonances for these high energy neutri- 
nos is therefore non-adiabatic, producing no significant v a — > v s conversion. The end result 
is a sterile neutrino energy spectrum as illustrated in Figure 1, with a characteristic mean 

rp jrp 

neutrino energy E approximately satisfying ^ f(E'/T)d(E'/T) « 8AL V J3 « 4C/3. As 
an example, for C = 0.02, we have E/T is 0.7, i.e., less than one-quarter of the average E/T 
for active species. The comoving free-streaming length of these sterile neutrinos becomes 

A,,.~40(^y7f^W (14) 



V30eV/ V 3 - 15 , 

which for a given sterile neutrino mass can be more than a factor of several shorter than 
that calculated from Eq. (Q). 

In this letter we limit our discussion of the sterile neutrino production to a temperature 
range below the Quark-Hardron phase transition, ~ 150 MeV. The reason for this is that 
the neutrino transformation process can be modeled confidently below this temperature. 
With this restriction, Eq. (|7|) shows that the pre-existing C has to be > 10~ 3 . The mass 
of the sterile neutrino Dark Matter candidate is then < 10 keV from Eq. (0). (C and the 
sterile neutrino mass are limited simultaneously because their product must yield Q u ~ 1.) 
If we also restrict the sterile neutrino mass to be m Vs > 10 2 eV in order to have potentially 
successful structure formation at galaxy scales, Eq. (|T3|) implies that £ < 10 _1 . For example, 
if the pre-existing 2L Ufi + L Ue +L Ut is 0.02, a sterile neutrino with m Va = 1160 eV and mixing 
with will yield Q u = 0.4 for h = 0.65. The comoving free-streaming length of this Dark 
Matter is only Af. s . ~ 0.4 Mpc, giving rise to a mass-scale cut-off of M ~ 1O 1O M , roughly 
the scale of dwarf galaxies. In this case, structure formation in CoolDM models differs from 
that in CDM models only at sub-galactic scales, such as in halo structure and formation 
history at high redshift. These differences may potentially be tested by observations of 



galaxy rotation curves [14 and high redshift galaxies [IJ],[L5fl. However, if the mass of the 



sterile neutrinos is as high as 10 keV, resulting from C ~ 10 -3 , then CoolDM essentially 
becomes CDM. It is intriguing that oscillation between active neutrinos and keV sterile 
neutrinos in supernovae in the presence of strong magnetic fields may potentially give rise 
to "pulsar kicks" [Tfj . 



An L between ~ 10~ 3 and ~ lO^ 1 necessarily implies that the asymmetries in one or 
more active neutrino species are of the same order (or larger if there is cancellation between 
asymmetries). Lepton asymmetries between ~ 1CT 3 and ~ 1CT 1 are not ruled out by Big 



Bang Nucleosynthesis |L7|]- In fact, if this asymmetry resides in the or v T sector, its 
impact on the primordial 4 He yield (which gives the constraint on lepton asymmetry) is 
well below the current detection level. If some of the asymmetry is in the v e sector, with a 
magnitude > 10~ 2 , its impact on the primordial 4 He yield is quite appreciable. Just where 
the limit on L„ e stands is hard to quantify precisely because of the large uncertainty in the 
primordial 4 He abundance measurements |jl|,[n|. In particular, Kohri et al. |2(J has argued 



that an L Ve « 0.015 is what is needed to bring a low set of primordial 4 He abundance values 
H| into agreement with the measured low primordial deuterium abundance [2T|. 



On the other hand, in individual neutrino sectors, lepton asymmetries of order ~ 10~ 3 
to ~ lO^ 1 are very large compared to the baryon asymmetry, ~ 10~ 10 . But this is not to say 
that their sum, the total lepton number asymmetry, cannot be as small as ~ 10~ 10 . In the 
absence of such a cancellation, however, the resultant large lepton number must be generated 
below the Electroweak phase transition to avoid the transfer between lepton number and 
baryon number. Models that generate a large lepton number despite a small baryon number 
asymmetry indeed exist [p^ , p3|| . 

In summary, we have proposed that ~ 10 2 eV to ~ 10 keV sterile neutrinos v s will be 
a Cool Dark Matter candidate if they were produced in the early universe via active-sterile 
neutrino resonant oscillation in the presence of a lepton number asymmetry. The amplitude 
of this lepton number asymmetry is ~ 10~ 3 to ~ 10 _1 for at least one active neutrino species, 
which is consistent with constraints from Big Bang Nucleosynthesis. The minimal set of free 
parameters in our Cool Dark Matter model, and its potentially testable predictions for halo 
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structures and proto-galaxy formation, all make the Cool Dark Matter scenario interesting 
and worthy of further invistigations. 

X. S. and G. M. F. are supported by NSF grant PHY98-00980 at UCSD. 
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Figure Captions: 

Figure 1. Illustration of the energy spectrum of the sterile neutrinos produced vs. the 
thermal spectrum of the active neutrino species. 
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